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ABSTRACT 

If an accretion disk around a black hole is illuminated by hard X-rays from non-thermal 
coronas, fluorescent iron lines will be emitted from the inner region of the accretion disk. The 
emission line profiles will show a variety of strong field effects, which may be used as a probe of 
the spin parameter of the black hole and the structure of the accretion disk. In this paper we 
generalize the previous relativistic line profile models by including both the black hole spinning 
effects and the non-axisymmetries of warped accretion disks. Our results show different features 
from the conventional calculations for either a flat disk around a Kerr black hole or a warped 
disk around a Schwarzschild black hole by presenting, at the same time, multiple peaks, rather 
long red tails and time variations of line proflles with the precession of the disk. We show disk 
images as seen by a distant observer, which are distorted by the strong gravity. Although we are 
primarily concerned with the iron K-shell lines in this paper, the calculation is general and is 
valid for any emission lines produced from a warped accretion disk around a black hole. 



Introduction 



Accretion power serves as an important source of energy in astrophysics, especially in galactic X-ray 
binaries (XRBs) and active galactic nuclei (AGNs). In these systems, black holes (BHs) are the possible 
candidates for the central accreting compact objects. Accretion onto both stellar-mass black holes in XRBs 
and supermassive black holes (10^ — IQ^Mq) in AGNs usually proceeds through an accretion disk, which 
may be truncated at the inner most stable circular orbit (ISCO) because of the general relativistic effect of 
black hole. Due to viscosity, plasma elements in the disk gradually spiral inwards to the black hole while 
losing their angular momenta until reaching the ISCO; they then plunge into the black hole, if there is no 
external force acting on them. 

The quasi-thermal spectra of the black ho le accretion disks a re relatively cold, which peak at optical/UV 
band in AGNs and soft X-ray band in XRBs (IFrank et alj|2002). If point- lik e hot coronae, above and below 
the center of the accretion disk (as in the lamp-post model (|Matt et al.lll99l[ )). illuminate the inner region of 
the disk by hard X-ray photons yielded via thermal Comptonization, iron atoms photoelectrically absorb the 
incident hard X-ray photons, then emit iron K-shell photons (centered at 6.40 — 6.97 keV, depending on the 
ionization state) through the fluoresc ent mechan i sm. I r on K-shell lin e s may be the strongest e mission lines 



i n the X -ray observatio n s of XRBs (IBarr et all (119851) . iMiller et all (|2002at ) for Cygnus X-lj^ 



Miller et al 



( 20Q2bl). iMiniutti et all (120041) . iRossi et al.l (120051) f or XTE J1650-500) a nd AGNs (T anaka et all (|l995[ ). 
Wilms et all |200l[ ). iBrenneman fc Reynolds! liooa ) . IMiniutti et Zl |2007t ) for the Seyfert-1 MCG-6-30-15, 
the most important source to date). Observations suggest that some two- horn broadened iron K-shell lines 
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originate from the innermost area of the accretion disk. The iron hnes could be a powerful tool to study the 
strong gravity within a few g ravitational radiu s = GM/c ^ and a probe of the black hole spin parameter 
(see reviews by iRevnoldsfc Nowak (,20031) andlMnid |2007l) V 



Fabian et al.l (|l989f ) modeled the relativistic line profile from a standard thin disk (se e lShakura fc Sunvaev 



(|l973l ) for Newtonian disk model: iNovikov &: Thornd (jl973) an d Page fc Thorne ( 1974 ) fo r gener al relativis- 
tic effects) orbiting around a Schwarzschild black hole. iLaoii (|l99l[ ) and I Bromley et al. J 1997) calculated 
the line profile from a standard thin disk orbiting around an extreme Kerr black hole. iFanton et al.l (jl997l ) 
calculated the line profile from a standard thin disk orbiting around a Kerr black hole with spin parameter 
a = cJ/GM (where J is the angular momentum and M is the mass of the black hole) as a free parameter, in 
their calculation the complication of computing azimuthal angle (j) has been avoided due to the axisymmetry 
of both the Kerr metric and the accretion disk. Thereafter, the line profiles from accretion disks with dif- 
ferent stru£turesJiaA^^ accretion disk orbiting around a Schwarzschild black 



hole (jHartnoU fc BlackmanlboOOt 



Wu et al 



20081). and for a flat accretion disk with spiral pattens orbi ting 



around a nonrotating black hole ([Hartnoll fc Blackmanll2002l ) or a rotating black hole (jKaras et al I2OO1I) . In 
these calculations, the accretion disks are optically thick and geometrically thin; furthermore it was assumed 
that the plas ma elements within ISCO contribute little to the line proflles, because they are h ot and almost 
fully i onized ( Young et al. 19981: Brenneman fc Reynolds 20061) . In contrast with the thin disk, Wu &: Wang 



( 2007 ) explored the role of disk self-shadowing effect due to significant geometrical thickness. iFuerst fc Wu 



20071 ) investigated the infiuence of geometrical thickness of relativistic accretion torus on the line profiles. 
It is interesting to m odel the relativistic ir o n K-s hell lines emerging from the warped disks orbiting 



around Kerr black holes. iBardeen fc PettersonI (jl975[ ) showed that for Kerr black holes the dragging of 



inertial frames can lead to coupling between the spin of the black hole and the orbital angular momentum 
of the disk, and that the component of the viscous torque, which is in the disk plane and perpendicular 
to the line of nodes, will tend to force the inner part of the disk to lie on the equatorial plane of the Kerr 
black hole. N evertheless, misalignment of the inner part of the di s k may take place even in the presence 



of dissipation. iDemianski fc Ivanovl (jl997[ ) and llvanov fc Illarionovl (|l997l ) showed that the misaligned disk 



adopts a steady warped shape in which the tilt angle is an oscillatory function of radius. iLubow et al.l (j2002l ) 
considered the time evolution of the warped disk under the condition that the warping propagates in a wave 
like manner, and concluded that, in low-viscosity disk around a Kerr black hole, the inner parts of the disk 



Pringle 


19961 


Malonev et al. 


1998) 



19961 ). In o bservations, evidence of warped disks in XRB includes the 35-days periodic variati on of X-ray flu x 
of Her X-1 ( Tananbaum et al. 1972 : Katz 19731) and the precessing relativistic jet in SS 433 (Margon 1984 ): 
for AGNs disk warping can e xplain the radial dependence of d eclination between high-velocity maser and 



systemic maser in NGC 4258 (jMiyoshi et al 



1995 



Mart: 



momentum of the radio jets and the host galactic disks ( Kinney et al. 200(]t Schmitt et al. 2002) 



2008f) and the misalignment between the angular 



The motivation of this work is multiple. Firstly, the iron K-shell lines carry the information about the 
spacetime around a rotating back hole; fitting line profiles from the observational data with better spectral 
resolution and signal-noise ratio by m ore realistic disk line mod els may provide more reliable estimations of 
the spin parameter of the black hole ([Abramowicz et al.l |2010|) , and distinguish the warping effects on line 
profiles from others. Secondly, the disk lines may be used to diagnose the geometry and kinematics of the 
disk itself, and provide observational signatures for the existence and characteristics (magnitude, twist-free 
or twisted) of the warping. 



The paper is organized as follows. In Sec. ID we review the properties of photon trajectories in the Kerr 
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spacetime, and summarize the relevant equations needed in our calculation. Next, we describe in Sec. |3ll 
the geometric formalism for the warped disk and in Sec. |3l2 the kinematic structure. These considerations 
are employed in Sec. SI where we show our main results, including the emission line profiles and the images 
of the warped disk around a Kerr black hole observed from a distant place. We summarize the main results 
and discuss their possible implications in Sec. [31 To allow the main ideas of the paper to be as clear as 
possible, several sets of details have been relegated to appendix A, in which we give the detailed procedure 
for deriving the expression of the ratio g for the wa rped disk. The relat ivistic calculations presented in this 
paper use the notational conventions of the text by Misner et al. ( 19731 ): in particular we use c = G = 1. 



2. Photon trajectories in the Kerr spacetime 

When calculating the relativistic line profiles, one needs to trace the trajectories of the photons in the 
Kerr spacetime. In the standard Boyer-Lindquist coordinates, the Kerr metric can be written as 

ds^ = -Y.AA-^df + sin^ eAY.-^ {d<j) - uidtf + EA^^dr^ + J^dO^ , (1) 

where 

S = + cos2 9 , (2) 

A = r^ - 2Mr , (3) 

A (r^ + a^) ^ - sin^ 6IA = S (r^ + a^) + 2AfaV sin^ 9 , (4) 

1-1 



uj^2aMrA-'. (5) 

This metric contains two parameters, the mass of the black hole M and the specific angular momentum 
or the spin parameter a — J /M . The greater root r+ = M + \/M'^ — of A = corresponds to the 
event horizon which is an one-way membrane in the Kerr spacetime. We notice that the right-hand-sides of 
Eqs.[2HS]are always positive outside the event horizon. The dragging of inertial frame forbids any observer 
remaining static within the static limit rsi{9) = M + ^/^P^^^^a?co^^ . 

Given a spacetime background, the motion of a photon obeys the null geodesic equation Vp • p = 0, 
where p is the 4-momentum of the photon and Vp is the directional derivative operator along p . This is 
a set of second-order ordinary differential equation, which can be solved in principle by the Runge-Kutta 
method with given initial position and momentum of the photon. Alternatively, in Kerr spacetime one can 
take advantage of its axisymmetric property (the corresponding constants of motion include the energy of 
the photon E, the angular momentum along the spin axis L^, and the norm of the 4-momentum 5), and solve 
the trajectory of photo n in the equat orial plane by the Lagrangian approach. For the general null geodesies 
in the Kerr spacetime, Carter ( 19681 ) demonstrated the separability of the Hamilton- Jacobi equation and 



discovered the existence the fourth constant of motion (i.e. Carter constant To obtain the equations 
governing the motion of the photons, one can take the partial derivatives of the Hamilton's principle function 



S^-5\-Et + L,(j) + / ^ ^ ' dr + / y/e(9)d9 , (6) 



with respect to E, L^, 5, and =S then set them to zero. After some manipulations (see Ichandrasekhai ( 19831 ) 



for details) , one can find that the equation governing the motion of the photon {6 = 0) in r and 9 is 

dr dO 
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where 



R{r, ^, r/) = + (a^ - - ??) + 2M (?? + (C - a)') r 
6(0, C, ?7) = ?7 + cos^ 9 - cot^ 6* , 



(8) 
(9) 



and P is the parameter along the photon trajectory. Here we have defined two dimensionless impact pa- 
rameters ^ = Lz/E and 77 = ^/E'^. The equations of motion for other coordinates (including the affine 
parameter A) can be expressed as combinations of the integrals about r and 6 



A 



r dr 



cos2 9d9 



t = E\ + 2M r (Er 



aLz 



a^E 



dr 



{{r^ + a^)E - aLz) 



dr 



A^Rir,C,v) 



A^R{r,^,r,) ' 
/ (Lz csc^ 9 - aE) 



de 



(10) 



(11) 



(12) 



The upper and lower limits in the integrals about r and 9 correspond to the turning points of the photon 
trajectories which are the roots of quartic Eqs. [8] and [9] with different choices of the impact parameters. 
The l imits are not written explicitly here due to the complexity (see Appendix A in iRauch fc Blandford 
(|l994i ) and Tables in lPexter fc AgoJ (|2009l ) for details). Eqs. H [H [Til and [Hare elliptical integrals, which 
ca n be reduced to th e Carlson's standard forms, then be evaluated numerically by the subroutines provided 
( 1992ri. In the ca lculation of photon trajectories, we use the eflScient subroutine geokerr 



m 

CT' 



Press et al. 



eated by iDexter fc Agoll (120091 ) which has integrated the necessary elements described above, to calculate 
all geodesic coordinates semi- analytically. In addition we shoot photons from the observer's plane of the sky 
backwards towards the central black hole and search for the intersections of the photons with the disk. This 
is a more efficient way than tracing the photons from the disk to the observer in our application. 



3. Geometric and kinematic structure of the warped accretion disk 

3.1. Geometric structure of the disk 

In our calculations, we use a series of concentric rings with increasing radii to describe the shape of the 
geometrically thin warped disk. The spherical (r, 0, </>) and cylindrical coordinate system {p, 0, h) are both 
appeared in our calculations, and they can be transformed from each other by 



and 



h — r sin /3 cos ip 
— p tan (3 cos tfj , 



1 /2 

7' (1 — sin^ /? cos^ (p) 
sin^/3 



(13) 



= r 1 - 



1 + tan^ -0 cos^ p 



(14) 



where we use p cos tp = r cos (p cos /3, tan ( p = tan V) tan and tjj = (t>—l in the derivation above. The formalism 
of the warped disk has been discussed by Terquem &: Bertoul ( 19931 ): HartnoU fc Blackman ( 2000 ): Wu et al 
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Fig. 1. — The geometry of the equatorial plane of the rotating black hole and the orbital plane of the plasma 
particle at r. The particle orbits around the central black hole along the circular ring in the orbital plane. 
The concentric rings with different /3 and 7 compose of the warped accretion disk. 



(|2008l ). Here we adopt the function expression of the two Eulerian angles used bv lWu et al.l (|2008l ). and the 
warping of disk is in an oscillation shape: 



J = lo + ni exp 7i2 



/3 ~ n3 sin 



Tout 



TT / r — 

2 V'^out 



(15) 



(16) 



We further specify that the inner radius of the disk r-m is always fixed at the ISCO, and the outer radius 
Tout at 50 rg. The parameters ni, n2, ns are used to describe the magnitude of warping. 70 measures the 
azimuthal viewing angle of the observer relative to the warped disk. Choosing different values for 70 is 
equivalent to changing the observational time during the disk precession. Therefore we can study the time 
variations of the line profiles and the disk images. 



3.2. Kinematic structure of the disk 

The kinematic structure of the warped disk carries the information about the motion of emitting particles 
in the vicinity of the Kerr black hole. It mainly contributes to the Doppler shifts (the only contributor to 
blue-shift) in the line profiles. We consider that each plasma element follows a quasi-Keplerian orbit in a 
concentric ring, which means the radial velocity is much less than the transverse velocity as in the standard 
thin disks, and the angular velocity of the spatially circular orbit can be written as 

n^{a + ey/r3/M)-\ (17) 

where e = 1 corresponds to the prograde orbit and e = — 1 the retrograde orbit with respect to the spinning 
of the black hole. This angular velocity is modified by the relativistic effect of the black hole by the inclusion 
of an additional spin parameter a. Strictly speaking, the circular orbit exists only in the equatorial plane. 
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However we assume that the tiled angle is small so that Eq. [17] is a good approximation for the warped 
disk. For the clarity in this paper, we simply list the 4-velocity of particles below and relegate the detailed 
derivation to Appendix A. 



^ fY,AA-^ - sin^ eAJ:-^ (j) - - T,eA , (18) 



2 \ -1/2 

SAA^i - sin^ OAJ:-^ {6-uj 

u^'^O, (19) 
— 9u^ = flu^ (— cos ip cos /3 sin ip cos 6 + sin ip cos (p cos 9 + sin /3 sin ip sin 9) , (20) 

u'^ — — (sin t/i cos /3 sin (/j + cos ^ cos v?) . (21) 

sm^ 



4. Emission line profiles and warped disk images 

4.1. Emission line profiles 

The emission line profile produced by an accretion disk around a black hole is influenced by Doppler 
shift, gravitational red-shift, beaming effect, and gravitational lensing effect. The changing of photon energy 
along a bundle of photon trajectories (geodesic congruence) can be characterized by the ratio g of the energy 
-Eobs measured by a local inertial (rest) observer at asymptotic infinity to the emitted energy iJom measured 
by an observer comoving with the plasma element on the accretion disk, 

5=|^ = ^^^^H^-(l + -r^ , (22) 

where z = (Aobs — -^cm)/''^cm is the red-shift as usually defined, u'^ in the denominator and numerator are 
the four- velocity of the emitting plasma element and the observer at infinity respectively, and is the four- 
momentum of the photon which propagates from the plasma element to the observer along a null geodesic. 
The total observed flux -Fobs is determined by integrating the specific flux dFobs over all the plasma elements 
on the accreting disk. For the line profile, the specific flux d-F'obs(-E'obs) at observed energy i?obs can be 
expressed as: 

dFobs(£'obs) = ^bs(-Eobs)df^obs , (23) 

where Jobs is the observed specific intensity and dfiobs is the observed solid angle subtended by each plasma 
element. Using Liouville's theorem 

/obs/l'is = lerJl^L (24) 

we find 

/obs /omi^obs/'^em = ^cmQ^ = /cm(l + z)'^ . (25) 

And the observed solid angle dfiobs is naturally the size of the pixels in the observer's screen. The detailed 
derivation of the factor g for the warped disk is shown in Appendix [Xj We ignore the intrinsic line width and 
assume Icm = £(5(i?cm — Eq), where Ep is the re st energy of photons, e ^ (r/M)'^ is the surface emissivity, 
and q is the emissivity index ( Fabian et al. IQSOf) . 



In the actual computation, we trace every null geodesic from the observer's screen backwards towards 
the central black hole and search for the interactions of the trajectories with the surface of the disk. The 
shadowing from the observer due to the warping of the disk is taken into account; we ignore the shadowing 
from the coronae, since in our calculation we assume that the scale height of the coronae is large compared to 
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the size of the disk of our interest, therefore this shadowing may not be important. We calculate the observed 
specific flux d-Fobs from each plasma element, and add its contribution into a energy bin according to the 
observed energy of photons (-Eobs — gEcm)- Then we can obtain a line profile by plotting count number 
(d-FbinZ-Ebin) vcrsus Fbin- The number of bins we choose determines the spectrum resolution. We assume an 
optically thick disk, thus we consider only the direct photons and neglect the photons which would circular 
around the black hole and strike the disk. 



4-. 1.1. Changing of the spin parameter a 

The overall properties of the line profiles can be understood by carefully examining the contribution 
of the specific flux from each of the narrow rings of the disk, they are determined by the spin parameter, 
the warping parameters, the emissivity index, shadowing and orientation of the disk relative to a distant 
observer (both the inclination angle and azimuthal angle). 

Figure [2] shows the line profiles from the warped accretion disk with the warping parameters ni = 0, 
rt2 = 1, and 713 = 0.95. The inner radius rin of the prograde disk (e = 1) locates at the innermost stable 
circular orbit (ISCO), which is 4.23 rg for the black hole spin parameter a — 0.5Af (maybe a good estimate 
for the supermassive black hole in Sgr A*). The outer radius is chosen as rout = 50 r^. From Eq. [TBI we 
know that the tilt angles are 0° and 54° for the narrow rings of the disk at the inner radius and the outer 
radius respectively. This warped disk is twist-free which means that the twist angle 7 does not depend on 
the radius [n.i = 0). The emissivity index q is taken to be —2. The horizontal axis is the g factor, i.e. the 
observed photon energy per unit iron K-shell photon energy, while the vertical axis gives flux in arbitrary 
units at different observed photon energies. In each panel, there are three colored line profiles representing 
the inclination angles of 10° (red), 50° (blue), and 70° (black) measured from the spin axis of the black 
hole respectively. Panels (a)-(h) contain the line profiles seen from different azimuthal viewing angles, for the 
retrograde precession of the warped disk about the spin axis of the black hole from 0° to 315° with uniform 
angular intervals of 45°. Panel (a) shows the = 0° case, which is defined as the observer viewing from the 
direction furthest from the lowest point of the disk. Panel (i) shows the line profiles from the standard flat 
disk with the same inner and outer radius for comparison. 

The line profile is sensitive to the inclination angle of the observer. For the low inclination angle = 10° 
(red), the hue profile in panels (a)-(h) tends to be a single broad peak with the same endpoints of the red 
and blue tails, and almost all of the photons are red-shifted due to strong gravitational red-shift and weak 
Doppler blue-shift. Those line profiles are similar to the fiat disk line profile in panel (i). With the increase 
of the inclination angle (blue and black), the two-peak or multi-peak structures start to emerge. Compared 
with the fiat disk case where the two peaks are shown (the blue peak is always brighter than the red one), 
the line profiles from the warped disk present considerable changing of the features, for example the red 
peak can be comparable or stronger than the blue one in some cases. The deviations are produced by the 
warping structure and shadowing effect, which can only be seen from the images of the disk in Fig. |3] (for 
6* = 50°). 

For the large inclination angles [0 = 50° and 9 = 70°), the shape of the line profiles is also sensitive 
to the azimuthal angle. This is because the observable regions of the disk are changing with the precession 
of the disk over a long period (see images in Figs. |3]for 6* = 50°). It appears that the line profiles are less 
sensitive to the low inclination angles {9 = 10°). Since at the low inclination, the inner region of the disk is 
almost face-on (for the twist-free warped disk the innermost area coincides with the equatorial plane) and 
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the shadowing effect due to the warping is less important. We note that the black lines in panel (a) {4> — 0°) 
and panel (e) {(/) = 180°) for 9 = 70° are similar, and the relatively narrower tails for (f) — 0° is due to the 
fact that part of the strongly blue-shifted and red-shifted area is shadowed. 

Figure |4] shows the line profiles from the accretion disk with the same warping parameters ni, n2, and 
77.3 as in Fig ure [H but for a more relativistic system, where the black hole has an extreme spin parameter 
a = 0.998M (jThorndll974 ). The inner radius ri„ for the accretion disk is only 1.23 Vg. The general features 



mentioned for Figure hold here. In addition, the effect of the gravitational red-shift and light bending is 
more significant. Some photon trajectories which are initially retrograde may be forced to turn back to be 
prograde by the frame dragging effect. In this more relativistic system, the red tails are rather long and 
extend to as far as g = 0.1, while the blue tails end at about the same places as in Figure El although they 
are cut off more sharply for largest inclinations {9 = 70°). However, the line profiles around g = 1 (the rest 
energy) are very similar to Figure [21 since the contribution to the flux here is mainly from the rings in the 
outer region of the disk (r > 10 r^), where the frame dragging effect is much smaller (w cx r~^). 

Figure |5l shows the line profiles from the accretion disk with the same warping and spin parameters 
as in Figures m except we choose the disk to rotate around the black hole in a retrograde orbit (e — — 1). 
The retrograde orbiting of the disk will increase its inner radius ri„ to a larger value 7.55 Vg, which reduces 
the contribution from the innermost part of the disk so that the red and blue tails become slimmer than in 
Fig. [21 In addition the patterns are anti-symmetric about 180° compared to Fig. [2l 



4-. 1.2. Changing of the emissivity index q 

The detailed irradiation law is unknown and the emissivity of the iron line in the inner region of the 
disk is poorly understood so far. For simplicity, we assume that the emissivity is a power-law function of 
radius with the emissivity index as the power-law index (i.e. e ^ {r /MY). 

Figure [6l shows the line profiles from the warped disk with the same warping parameters and the spin 
parameter as in Figure [21 The difference is that the inclination angle 9 in all cases here are chosen as 70° 
measured from the black hole spin axis. In each panel, the three colored lines represent the line profiles with 
the emissivity index q being —2 (red), —2.5 (blue), and —3 (black) respectively. Panels (a)-(h) contain the 
line profiles seen from different azimuthal viewing angles as in Figs.[2H5l and panel (i) shows the line profiles 
from the standard flat disk for comparison. All the line profiles are normalized to unity. 

Generally speaking, smaller emissivity index indicates relatively larger contribution from the inner region 
of the disk, if we do not consider the shadowing effect which prevents part of the inner region to be seen 
in some cases. The inner regions contribute mostly to the flux at blue and red tails while the outer region 
to that around 5 = 1 (the rest energy). Therefore, for smaller emissivity index (black than red), the flux 
in the red/blue tails is more intensive than that for larger emissivity index. However the shadowing effect 
is pronounced in some cases. For example in panel (a) (0 = 0°), part of the most blue-shifted region is 
unobscrvable due to the shadowing, the blue tails for three different q are almost identical. Moreover, they 
are less blue-shifted [g ~ 1.2) and are not steep as others. Similar feature appears in panel (b) ((/> — 45°) 
and panel (h) [cf) — 315°) too. In panel (e) ((/) = 180°), the disk is fully observed (least shadowing), the main 
feature of the line profile is similar to the flat disk where there is no shadowing at all. In panels (c)-(g), 
where the shadowing effect is less important, the blue tails for smaller emissivity index are steeper, and it 
is a common feature shared by both the warped disk and the flat disk. Since the nature of the central black 
hole, the dynamic structure of the disk and the orientation to the observer are the same for profiles in each 
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panel, the blue and red tails representing different emissivity index end at the same places. 

Figure [7] shows the line profiles from the warped disk with the same warping parameter and the spin 
parameter as in Figure [Bl but the inclination angle 9 here are chosen as 50°. The features mentioned for 
Figure |6] remain true here. Furthermore the blue tails are reduced to smaller values {g < 1.2) and cut off 
more steeply due to smaller Doppler blue-shift for this smaller inclination angle; while the red tails end 
at about the same value {g ~ 0.4), because they are mainly from the gravitational red-shift which is not 
sensitive to the inclination angle. 

Figure [8] shows the line profiles from the warped disk with the same warping parameter and inclination 
angle as in Figure[7]but with the spin parameter a = 0.998Af. For such extreme spin rate, the inner radius of 
the accretion disk can be very close to the event horizon (r^, — 1.0632 r^), where the gravitational red-shift 
dominates. Therefore the red tails in Figure |8] are much longer than those in Figure [T] For the reason we 
motioned above, smaller emissivity index usually produces more intensive red tails, the red tails for q = —3 
(black) is stronger than those of the other two lines. In extreme case, for example panel (a), a bump appears 
in the red tail oi q = —3. 



Discussion and Conclusions 



To date there are two practical methods for measuring the spin of astrophysical black holes, namely 
the X-ray continuum fitting and the relativistic iron lines. T he continuum-fitting method so far can only be 
applied to stellar- mass black holes (|McClintock et al.ll201l[ ). whereas the relativistic iron lines can be used 
as a probe to both the spins of stellar-mass black holes in X-ray binaries and supermassive black holes in 
active galactic nuclei. Furthermore if the accretion disks in these systems are warped, it is possible that the 
relativistic iron line profiles may serve as a way to diagnose the structures of the disks. 

In this paper we have developed a method for calculating the emission line profiles from the warped 
accretion disks orbiting around rotating black holes. It essentially generalizes previous calculations by 
including the black hole spin, the disk warping and the shadowing effect at the same time. We have also 
presented the disk images which can help us understand the general features of the line profiles from different 
disk and black hole systems. The detailed calculation concerns the computation of the trajectories of the 
photons emitted from the disk and the velocity field of the plasma elements orbiting around the center 
black hole. The formalism described above makes no assumption about the emitting energy of the photon. 
Therefore it can be applied to any emission lines which could be detected in the immediate vicinity of the 
black hole, such as the carbon, nitrogen, and oxygen lines although they are not as prominent as the iron 
lines. 

We assumed a set of warping parameters for the geometry of the twist-free warping disk and Keplerian 
orbit for the motion of the plasma elements as examples. Even for this simplified system we have found a 
number of interesting phenomena: (1) Different from the two-peak feature generally found in flat disks, the 
line proflles from the warped disk show multiple peaks. (2) The line profiles for some cases present sharper 
red tails and/or s ofter blue tails, whi ch may provide an alter native explanation to observations for some 
Seyfert 1 galaxies (jNandra et al.lll997l ) and Seyfert 2 galaxies (jTurner et al.lll997l ). (3) When the disks are 
orbiting around a highly spinning black hole, a rather long red tail is a common feature shared between flat 
disks and warped disks, if the shadowing effect is not important. (4) At low inclinations, the line proflles from 
different azimuthal angles differ very little from each other, indicating that the warping has little influence 
on the line proflle and it may not provide much information for diagnosing the disk structure. (5) The line 
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profile is also sensitive to the illumination law: smaller emissivity index q generally indicates stronger red 
and blue tails. In some extreme cases, red bumps may grow in the red tails. (6) It may not be true that 
a sharp blue tail is a common feature in all axisymmetric disk m odels, and is expected for time-integrated 
profiles even if the disk has strong asymmetry or inhomogeneities ( Bromlev et al. IQQSl) . The shadowing we 
encounter here may cause a soft blue tail in the integrated profile. (7) With the precession of the warped 
disk, time variations of line profiles may be present, and may be a possible observational signature for the 
warped disk. (8) Photometric variations of the line fluxes may be due to the varying amount of area of the 
disk facing to the observer. 

It is nontrivial to constrain the model parameters by fitting the theoretical line profiles to observations. 
However, from the above results and discussions we can see that generally the red tails of the line profiles 
provide the information about the spin of the black hole, while the blue tails the inclination angle. And 
the warping would cause complicated multiple peaks (which depends on the warping parameters), and time- 
variation of line profiles with the precession of disk. In the large inclination cases, the warping would cause 
significant disk shadowing which may prevent all/part of the inner area to be seen. Smaller emissivity index 
q normally indicates more intensive red and blue tails in the line profiles. Our present work is an initiative 
of a systematic investigation about the impact of the disk warping and shadowi ng on the line profile. It can 
be considered a s an extension of previous work on line profiles for warped disk ( Hartnoll fc Blackman 20001 : 
Wu et al. 20081) . Our line profiles distinguish from the line profiles from non-rotating black hole systems 
by presenting longer red tails which could extend to g — 0.1 for highly spinning cases. And the overall 
line feature is smooth and has simi l arity to the flat disk if the shadowing is not significant. These are 



different from iHartnoll fc BlackmanI (|2000l ). in which there are abrupt spike-like structures in some cases. 
We provide an altern ative explanation to the variation of line profile which is complementary to the picture 



of lKaras et al.l (|200lh . who investigate the non- axisymmetric patterns in the disk surface. However the non- 
axisymmet ric patterns are un likely to produce mul tiple peak structure . Due to the axisymmetry, both the 
thick disk (jWu fc Wanall2007l ) and accretion torus (jFuerst fc Wull2007l ) would not show the time- variant of 
line profile with the precession and the line profiles are usually single-peaked and double-peaked. 

Long-period photometric variation of X-ray is observed for many accreting systems. The precession 
of warped disk may be an explanation for this pheno menon. The time-scale characterizing the precession 
induced by tidal force of the companion star in XRB ( Wiiers fc Pringl3 19991 ) is 



Porh 



M 



1/2 



R 



101 



-3/2 



(26) 



where Porb is the orbital period of the binary, and q is the ratio of the mass of the companion star to the 
central accreting obj ect. The time-scale characterizing the precession induced by radiation can be found in 
Malonev et al. ( 19961) . The typical precession period for an AGN may be much longer than that of the XRB. 
The purpose of showing this formula is to point out that apart from the time-averaged line profile (with 
typical explosion time longer than 10^ s), it is possible to study the time variation of the iron line profiles 
on the precession timescale with the next generation space X-ray observatories which have sufficiently large 
collecting area. 

With the advance of the millimeter /sub- millimet er very-long baseline i nterferometer (VLBI), the pro- 
posed Earth-based Event Horizon Telescope (EHT) ( Doeleman et al. 20091) and the on-going space- VLBI 
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(such as vsofQ and RadioAstroro) promise to provide microarcsecond imaging resolution, which is suffi- 
cient to resolve the event horizons and the inner regions of the accretion disks around a handful of super- 
massive black holes. The calculated disk images may help understand the possible asymmetric intensity 
profile, blue/red-shift, and other observational signatures caused by the warping of the accretion disk and 
the rotating of the black hole. In the following work we will try to constrain the model parameters by fitting 
the theoretical profiles of the predicted spectra to actual data. 

Although we have given results for only one set of warping parameters, our formalism can be easily 
extended to more general disks. Work is underway on investigating systems with different warping parame- 
ters. By using the method of moments, we can explore the large warping parameter space (ni, rt2, n^) and 
study the distribution of the warping parameters in the moment space. In practice, it may be helpful to 
compare the positions of the prominent line peaks in the warped disks with the corre spondin g peaks in the 



flat d isks. And it would be easier than fitting the overall line profiles (jMurphv et al . 2009; .Sochora et al 



20111 ). Through this approaches we may connect the calculated line profiles with the current and/or future 



observations. 
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A. Expression of g for the warped disks 

This section describes how to derive the expression of g, the ratio of observed photon energy to the 
emitted energy for the warped disk. 



A.l. 4- velocity of the particle 

In each concentric ring, we let the angular position of the each plasma element he (p — fit, with t the 
coordinate time in the Kerr metric and the plasma particle's velocity components are 

v' = dxi' /dt — V {— sin (p, cos (p, 0) , (Al) 

where v — rfl — r{a + eyVYAf)^^, e = 1 is for co-rotating/prograde orbiting, e =^ —1 is for counter- 
rotating/retrograde orbiting. In addition, we need to transfer the velocity into the coordinate frame in 
which the z axis is aligned with the black hole spin. In Figure [TJ the inclined x'y' plane is the orbital plane, 
and the y' axis coincides with the line of nodes and points at the ascending node. The xy plane is the 
equatorial plane which has a tile angle /3 relative to the x'y' plane. The projection of the z' axis on the xy 



^http:/ /www. vsop.isas.ac.jp/ 

^http: / /www. asc.rssi.ru/radioastron/ 
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plane has a twist angle 7 relative to the x axis. We find the components Vi in the xyz coordinate basis using 
Vi = TijiVji, where 

(cos 7 cos (3 — sin 7 — cos 7 sin (3 \ 
sin 7 cos /3 cos 7 — sin A sin /3 , (A2) 
sin/3 cos/3 / 

from which we get 

Vx = V {— COS 7 cos P sin ip — sin 7 cos (^) 

Uj^ = w (— sin 7 cos /3 sin 1^9 + cos 7 cos (ys) (A3) 
= w (— sin /3 sin Lp) . 

We then calculate the velocity components in the spherical coordinate by projecting the velocity v onto the 
spherical basis e^, eg, 

Vr = Q 

vg = V (— cos ip cos (3 sin cp cos 9 + sin t/j cos p cos + sin /3 sin c/j sin 9) (A4) 
= V (sin V' cos /3 sin (/? + cos V' cos (^3) . 

where ip — (j) — j is the azimuthal angle of the plasma element measured from the projection of the x' axis 
on the xy plane. 

In the Kerr spacetime, if we set r as the proper time of the plasma element, it satisfies 

- dr^ = ds^ = -SAA-^t^ + sin^ (dcf) - ujdtf + SA^^dr^ + ^9"^ , (A5) 

and we divide both sides by dt^ 

dr^ / ■ \ 2 

- — = -Y.AA~^+&m^9AY.^^[(i,~ujj +T.9^ . (A6) 

Therefore we can write the 4- velocity of the plasma element as follow 



2 \ -1/2 

sin^ME-i ((i-cjl 



dt 

d^ ^ 


(^EAA 


dr 


dr dt 


d^ ^ 


dt dr 


d9 


d9 dt 


d^ ^ 
dcp 


dt dr 
d(f> dt 


d7 ~ 


dt dr 



dt dr 

u- = — = — — "-"^^ 

r 

^ d(p d<p dt -0 Q 
u — -r- — — :-— = (pu — u — r—„ 

r sm 9 



(A7) 



which are identical to Eq. [18] in Sec. 13.21 if insert Eq. IA4I into Eq. IA7I We note that the radial component 
of the velocity is zero because of we have assumed that the particle follows Keplerian orbit. And the term 
under the square root of Eg. IA7l is always positive if r > ngco- 



A. 2. expression for g 

The 4- velocity of the particle are expressed by the Boyer-Lindquist coordinates in Eq. IA7[ we then 
derive the expression for g. As defined 

^^1^^ (A8) 
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The observer rest at infinity has 4-velocity p^=(-l, 0, 0, 0), thus in the numerator (u'*p^)obs 
the denominator 

(w''P/.)em = -Eu° ± Eu°Veve/r + L^u°v^/ sine , 

therefore we get 

1 

^ ^ uO{-l±V&ve/r + ^v^/sme) ' 
where the sign is determined by the zenithal emitting direction. 
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Fig. 2. — Line profiles from a warped disk with ni = 0, n2 = 1, na = 0.95, ri„ = 4.23 rg and rout = 50 rg. 
The spin parameter a of the Kerr black hole is 0.5M. The emissivity index q is taken to be —2. The 
horizontal axis is the g factor, i.e. the observed photon energy per unit iron K-shell photon energy; the 
vertical axis is flux in arbitrary units. In each panel, three colored line profiles represent inclination angles 
of 10° (red), 50° (blue), and 70° (black) measured from the spin axis of the black hole respectively. Panels 
(a)-(h) contain the line profiles seen from different aziniuthal angles cj) of 0°, 45°, 90°, 135°, 180°, 225°, 
270° and 315°. Panel (a) shows the = 0° case which is defined as the observer viewing from the direction 
furthest from the lowest point of the disk. Panel (i) shows the line profile from a standard fiat disk for 
comparison. 
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Fig. 3. — Images of a geometrically thin and optically thick warped disk around a Kerr black hole. The 
spin parameter a is 0.5M, and the observer's inclination angle Q is 50° measured from the spin axis. The 
warping parameters are the same as previous figures (n\ = 0, n2 = 1, "-3 = 0.95). The inner and outer 
radius are 4.23 Tg and 50 Vg. a is the horizontal axis of the observer's photographic plate, and /3 is the 
vertical axis. The false color contour maps show the ratio of the observed energy to the emitted energy. The 
blue shaded areas represent the regions of the disk where photons emitted are blue-shifted, while the red 
shaded areas represent the regions where photons are red-shifted. The large areas colored with cyan and 
green represent the regions where photon energies are weakly blue and red-shifted. The white areas are the 
zero-shift regions, where the gravitational red-shift is balanced with the Doppler blue-shift. Panels (a)-(h) 
contain the line profiles seen from different azimuthal angles (j) of 0°, 45°, 90°, 135°, 180°, 225°, 270° and 
315°. Panel (a) shows the = 0° case which is defined as the observer viewing from the direction furthest 
from the lowest point of the disk. Panel (i) shows the disk image from a standard flat disk for comparison. 
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Fig. 6. — Line profiles from a warped disk with, ni = 0, 712 = 1, = 0.95, ri„ = 4.23 and rout = 50 r^. 
The spin parameter a of the Kerr black hole is 0.5M. The emissivity index q is taken to be —2. The 
horizontal axis is the g factor, i.e. the observed photon energy per unit iron K-shell photon energy; the 
vertical axis is flux in arbitrary units. The inclination angle 6 is 70° measured from the spin axis of the 
black hole. In each panel, three line profiles represent emissivity index q of —2 (red), —2.5 (blue), and —3 
(black) respectively. Panels (a)-(h) contain the line profiles seen from different azimuthal angles 4> of 0°, 
45°, 90°, 135°, 180°, 225°, 270° and 315°. Panel (a) shows the = 0° case which is defined as the observer 
viewing from the direction furthest from the lowest point of the disk. Panel (i) shows the line profiles from 
a standard flat disk for comparison. 
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